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Using numerical relativity as guidance and the natural flexibility of the effective-one- body (EOB) 
model, we extend the latter so that it can successfully match the numerical relativity waveforms of 
non-spinning binary black holes during the last stages of inspiral, merger and ringdown. Here, by 
successfully, we mean with phase differences 8% of a gravitational- wave cycle accumulated until 
the end of the ringdown phase. We obtain this result by simply adding a 4 post-Newtonian order 
correction in the EOB radial potential and determining the (constant) coefficient by imposing high- 
matching performances with numerical waveforms of mass ratios mi/ m 2 = 1,2/3, 1/2 and = 1/4, 
mi and m 2 being the individual black- hole masses. The final black-hole mass and spin predicted 
by the numerical simulations are used to determine the ringdown frequency and decay time of three 
quasi-normal- mode damped sinusoids that are attached to the EOB inspiral- (plunge) waveform at 
the light-ring. The accurate EOB waveforms may be employed for coherent searchs of gravitational 
waves emitted by non-spinning coalescing binary black holes with ground-based laser- interferometer 
detectors. 

PACS numbers: 


I. INTRODUCTION 

The network of ground-based laser-interferometer 
gravitational- wave (GW) detectors, such as LIGO [1], 
VIRGO [4], GEO [2] and TAMA [3], are currently op- 
erating at design sensitivity (except for VIRGO which 
plans to reach design sensitivity within one year) and 
are searching for GWs in the frequency range of 10-10 3 
Hz. Within the next decade those detectors will likely 
be complemented by the laser-interferometer space an- 
tenna (LISA) [5], a joint venture between NASA and 
ESA, which will search for GWs in the frequency range 
5 x 10" 5 -10“ 2 Hz. 

Binary systems composed of black holes (BHs) and/or 
neutron stars (NSs) are among the most promising 
sources. The search for GWs from coalescing binary 
systems and the extraction of parameters are based on 
the matched-filtering technique [6] , which requires accu- 
rate knowledge of the waveform of the incoming signal. 
Recent comparisons [8-11] between numerical and post- 
Newtonian (PN) analytic waveforms emitted by non- 
spinning binary BH systems suggest that it should be 
possible to design purely analytic templates with the full 
numerics used to guide the patching together of the in- 
spiral and ringdown waveforms. This is an important 
avenue to template construction as eventually thousands 
of waveform templates may be needed to extract the GW 
signal from the. noise, an impossible demand for numeri- 
cal relativity (NR) alone. 

The best-developed analytic method for describing the 


two-body dynamics of comparable-mass BHs and pre- 
dicting the GW signal is undoubtedly the PN method [7], 
which for compact bodies is essentially an expansion in 
the characteristic orbital velocity v/c. Predictions are 
currently available through 3.5PN order [14-17] (t? 7 /c 7 ), 
if the compact objects do not carry spin, and 2.5PN or- 
der [18] ( v 5 /c 5 ) if they carry spin. Resummation of the 
PN expansion aimed at pushing analytic calculations un- 
til the final stage of evolution, including the transition 
inspiral-merger-ringdown, have been proposed. In 1999, 
the authors of Ref. [21] introduced a non-perturbative re- 
summation of the two-body conservative dynamics, the 
so-called effective-one-body (EOB) approach [21]. The 
original EOB model was computed using the 2PN conser- 
vative dynamics. It was then extended to 3PN order [23] 
when the 3PN calculation was completed [17] and then 
to spinning BHs [24]. The EOB approach has been the 
only analytic approach able to predict, within ~ 10% of 
accuracy, the spin and mass of the final BH [25], and 
to provide a complete waveform, from inspiral to ring- 
down, for non-spinning [25, 26] and spinning, precessing 
binary systems [27]. To include accurately the radiation- 
reaction contribution, the EOB approach uses the Pade 
resummation of the GW flux, as proposed in Ref. [28]. 

As we shall see below, because of the reduction of the 
dynamics to a few crucial functions determining the in- 
spiral evolution, and because of the rather simple pro- 
cedure to match the inspiral (-plunge) waveform to the 
ringdown waveform, the EOB model is particularly suit- 
able for fitting to the numerical results [12, 13]. In this 
paper we shall employ its flexibility to obtain accurate 


2 


waveforms for non-spinning binary BHs with mass ratios 
ranging between mi/ m 2 = 1 and mi/ m 2 = 1/4, with mi 
and m 2 being the BH masses. The method also allows us 
to predict the waveforms for mass ratios mi/ m 2 <1/4. 
Those waveforms will be tested against numerical results 
when accurate long numerical simulations for mass ra- 
tios mi/ m 2 < 1/4 become available. In this paper the 
comparisons are carried out using simulations from the 
NASA-Goddard group. 

The paper is organized as follows. In Sec. II we briefly 
review the EOB model. In Sec. Ill we improve the EOB 
model by adding a 4PN order unknown coefficient to the 
two-body conservative dynamics. In Sec. IV we compare 
this improved EOB model to two accurate, long numer- 
ical simulations having mass ratios 1 : 1 and 1:4, de- 
termine the best-fit. 4PN order coefficient and discuss the 
matching performances. Section V summarizes our main 
conclusions. In Appendix A we extend the comparison to 
several multipole moments, and in Appendix B to shorter 
numerical simulations having mass ratios 1 : 2 and 2:3. 

II. THE EFFECTIVE- ONE-BODY MODEL 

The resummation technique discussed in this section, 
the EOB approach [21], was originally inspired by a sim- 


ilar approach introduced by Brezin, Itzykson and Zinn- 
Justin [22] to study two electromagnetically interacting 
particles having comparable masses. The basic idea of 
the EOB approach is to map the real conservative two- 
body dynamics up to the highest PN order available, 
onto an effective one-body problem, where a test par- 
ticle of mass (i = mi?7i2/M, with mi, m 2 the black-hole 
masses and M = mi + m 2 , moves in some effective back- 
ground metric g This mapping has been worked out 
within the Hamilton- Jacobi formalism, by imposing that 
whereas the action variables of the real and effective de- 
scription coincide, i.e. L rea i = Xreal = T e ff, where 
L denotes the total angular momentum, and X the ra- 
dial action variable, the energy axis is allowed to change, 
Ereai = f(E e ff ) , where / is a generic function determined 
by the mapping. By applying the above rules defining 
the mapping, it was found [21] that as long as radiation- 
reaction effects are not taken into account, the effective 
metric is just a deformation of the Schwarzschild metric, 
with deformation parameter 77 = p/M. 


The explicit expression of the non-spinning EOB effec- 
tive Hamiltonian through 3PN order is [21, 23]: 


#eff (r , p) = p £T e ff (r , p) = p A A(r) 


l + P 2 + ~ ( n ’P) 2 + ~2 Oi(P 2 ) 2 + 22 P 2 (n-p) 2 +z 3 (n-p) 4 ) 


( 1 ) 


with r and p being the reduced dimensionless variables; motion is constrained to a plane. Introducing polar co- 
il = r/r where we set r = |r|. In absence of spins the ordinates (r, <p,p r ,jp^), the EOB effective metric reads 

I 


dslfi = g e Jl dx M dx v = -A(r) c 2 dt 2 + dr 2 + r 2 (d9 2 + sin 2 9 dp 2 ) 


( 2 ) 


The EOB real Hamiltonian is 


arbitrary, but subject to the constraint 


ffreal - M H rea l = M ^ 1 + 2t) ^ — M , (3) 

Remarkably, as originally observed in Ref. [21], the 
mapping between the real and the effective Hamiltoni- 
ans given by Eq. (3) coincides with the mapping ob- 
tained in the context of quantum electrodynamics in 
Ref. [22], where the authors mapped the one-body rel- 
ativistic B aimer formula onto the two-body energy for- 
mula. Moreover, Eq. (3) holds at 2PN and 3PN or- 
der [23]. The coefficients 21,22 and £3 in Eq. (1) are 


821 4- 4 z 2 + 323 = 6(4 - 3t]) rj . (4) 


The coefficients A(r) and D(r) in Eq. ( 1 ) have been calcu- 
lated through 3PN order [ 21 , 23]. In the Taylor-expanded 
form they read: 


4 PN W = i-;+^+ 




rj- zi 


1 

r*’ 


(5) 


D ¥ N ( r ) ~ 1 - + \J z i + 2:2 + 277 (377 - 26)] L . (6) 
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FIG. 1: In the left panel we show the position of the ISCO and light-ring as function of the parameter w , for different mass 
ratios. In the right panel we show the energy for circular orbits as a function of the frequency evaluated from the EOB 
Hamiltonian and the radial potential as a function of the radial coordinate for a massless particle in the EOB model. The 
various curves refer to different PN orders. 


In principle we could explore the possibility of determin- 
ing some of the Zi coefficients through a fit with the nu- 
merical results. However, here we do not follow this pos- 
sibility and, as in previous works, set z\ = Z2 =0, thus 
z 3 = 2(4 — 377)77. The EOB effective potential A|? N (r) 
does not lead to a last-stable circular orbit (LSO), con- 
trary to what happens in the 2PN-accurate case [21]. 
This is due to the rather large value of the 3PN coeffi- 
cient 94/3 — 41/327T 2 ~ 18.688 entering the PN expansion 


of A(r). Replacing the PN-expanded form of A(r) by a 
Pade approximant cures this problem [23]. The Pade 
approximant is 


i 2PN 


at 2PN order and 


(r) 


r (—4 + 2r + 77) 
2r 2 + 277 + r 77 


(7) 


.3PN/ s r 2 [(0.4(77, 0) + 877 - 16) + r (8 - 2 r?)] , g . 

p s ( - 7 ' r 3 (8 - 2 rf) + r 2 [a 4 (??, 0) + 477] + r [2 a 4 (ry, 0) + 877] + 4[ry 2 + a 4 (r?, 0)] ’ v 


at 3PN order where 


and t = t/M, Q = u M [25], as [54] 


a 4 (?7, zi) 


94 

3 


41 

32 1 


zi 


(9) 


dr 

dt 

dip 

dt 

dp r 

dt 

dpy 

dt 


dH (rn n \ 

d P A Pr ' Plp) ' 

(10) 

~ oh. , 

w = -T (r,Pr,P v ) , 

OPv 

(11) 


(12) 

F v [u(r,p r ,p v )}, 

(13) 


To include radiation-reaction effects we write the EOB 

Hamilton equations in terms of the reduced quantities H where for the p component of the radiation-reaction force 
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FIG. 2 : Equal- mass binary. In the left panel we plot the NR and EOB frequencies and amplitudes, and the phase difference 
between the EOB and NR /122. In the right panel we compare the EOB and NR hii extracted perpendicularly to the orbital 
plane. We maximize only on the initial phase and time of arrival. 


we shall the P-approximant [25, 28] 


III. IMPROVING THE EFFECTIVE-ONE-BODY 
MODEL 


32 7 fPsi^V)' 

5 V u 1 - v u /v po \ e {ri) ’ 
(14) 

where v w = Q 1 ^ = ( d<p/dt ) 1 / 3 . The coefficients fp N 
can be read from Eqs. (XXX) in Ref. [42], while we set 
'Up oie = 1. In Refs. [44], the authors pointed out that a 
more accurate expression of the radiation-reaction force 
should not use the Keplerian relation between r and uj 
when the binary evolves inside the LSO. We find that 
this modification has a negligible effect on the waveform. 

The last crucial ingredient of the EOB model is the in- 
clusion of the ringdown phase. After the two BHs merge, 
the system settles down to a Kerr BH and emits quasi- 
normal modes (QNMs) [49, 50]. In the test-mass limit, 
77 <C 1, Refs. [50, 51] realized that when a test parti- 
cle falls radially below 3 M (the unstable light-ring of 
Schwarzschild), it immediately triggers the production 
of QNMs. In the equal-mass case rj = 1/4, Ref. [52] 
proposed the so-called close-limit approximation, which 
consists in switching from the two-body description to 
the one-body description (perturbed-BH) close to the 
light-ring location. Based on these observations, Ref. [21] 
modeled the merger as a very short (instantaneous) phase 
and matched the inspiral (-plunge) waveform at the light- 
ring to a damped sinusoid. The frequency and decay 
time were computed estimating the final BH mass and 
spin from the EOB energy and angular momentum at 
the matching point. The matching procedure has then 
been improved, by adding more QNMs, extending it to 
several multipole moments [8, 53], and applying it in an 
time-interval instead of one point in time [46] . 






Previous investigations [8, 10] focusing on compa- 
rable mass binaries, have shown that a non-negligible 
dephasing can accumulate at the transition inspiral(- 
plunge) to ringdown between the EOB waveform, com- 
puted through 3.5 PN order, and the NR waveform. The 
dephasing is caused by the much fast increase of the GW 
frequency in the 3.5PN-EOB model than in the NR simu- 
lation. Although the dephasing will prevent to accurately 
determine the binary parameters, it will not prevent to 
detect the waveform. In Ref. [10] the EOB matching 
to ringdown employed three extra parameters describing 
the time of matching, and the difference between the final 
BH mass (spin) and the energy (angular-momentum) at 
the point of matching. In this paper we shall improve the 
matching getting rid of the three extra parameters and 
reducing the dephasing to <0.08 of a GW cycles, thus 
providing accurate waveforms depending only on the BH 
masses to be used for coherent detection and parameter 
estimation. 


To decrease the differences between the EOB and NR 
waveforms during the last stages of inspiral and plunge, 
we introduce a 4PN order term in the effective potential 
A(r) and Pade approximate it using the approximant 
A P . 1. Similar modifications were employed in Ref. [33] 
to obtain better matches of the EOB model to quasi- 
equilibrium initial-data configurations [34]. Since the 
4PN term has not been calculated in PN theory, we shall 
denote it as “p4PN”, where “p” stands for pseudo . We 
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FIG. 3 : Binary with mass ratio 1 : 4 . In the left panel we plot the NR and EOB frequencies and amplitudes, and the phase 
difference between the EOB and NR /i22- In the right panel we compare the EOB and NR /122 extracted perpendicularly to 
the orbital plane. We maximize only on the initial phase and time of arrival. 


have 



Num( J 4pt PN ) 
Den(^ PN ) ’ 


with 


(15) 

I 


Num(Ap 4PN ) = r 3 [32 — 2 Arj — 404(77, 0) — 05(77,70)] + t -4 [04(77, 0) — 16 + 877] , (16) 

Den(j4p 4PN ) = -<14(77,0) - 805(77,70) - 804(77,0)77 + 205(77,70)77 - I677 2 + r [-804(77,0) - 405(77,70) - 204(77,0)77 - 

r A 

I677 2 ] + r 2 [-404(77,0) - 205(77,70) - I677] + r 3 [-204(77,0) - 05(77,70) — 877] + r 4 (-16 + 04(77,0) +877) , 

(17) 


where 

a 5 (r),w) =wt], (18) 

and w will be determined by the comparison with nu- 
merical results. In the left panel of Fig. 1 we show how 
the pseudo 4PN order term w modifies the position of 
the EOB LSO and light-ring (the last unstable orbit for 
a massless particle) for several binary mass ratios. Later 
on we shall see that the value of w that best fits the NR 
results (see Sec. IV) is w = 60. It always guarantees the 
presence of a LSO and a light ring. In the right panel 
of Fig 1, we show the circular-orbit energy computed 
with the EOB Hamiltonian, and the radial potential for 
a massless particle, at different PN orders having fixed 
w = 60. We notice that the LSO energy = 

—0.0185) and frequency (Mcj^pn = 0.1047) at p4PN or- 
der are closer to the corresponding values obtained using 
the 3PN-Taylor-expanded model for quasi-circular adia- 


batic orbits [29], and to the quasi-equilibrium initial-data 
approach [30] (see Fig. 16 and Table II in Ref. [30]). 
This could be a pure accident. In fact, it should be 
kept in mind that the LSO frequency computed from 
the 3PN-Taylor-expanded conservative dynamics is [29] 
McjgpN ~ 0.129 (Ej PN /M = —0.0193), quite close to 
the formation of the common apparent horizon in the 
NR simulation, and quite far from the frequency ~ 0.08 
at which the blurred plunge occurs [8]. What we can 
certainly say is that the EOB conservative dynamics at 
p4PN order is closer than at 3PN order to the 3PN 
Taylor-expanded conservative dynamics of quasi-circular 
adiabatic orbits [29]. 
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FIG. 4: For different mass ratios, we show how the FF (com- 
puted using white noise) depends on the parameter w. The 
vertical line refers to the value w = 60 which we employ in 
all subsequent analyses. 


IV. ACCURATE EFFECTIVE-ONE-BODY 
WAVEFORMS FOR INSPIRAL, MERGER AND 
RINGDOWN 

We now compare the p4PN EOB model of Sec. Ill 
with a few numerical simulations produced by the NASA- 
Goddard group. The numerical simulations compute 
directly the Weyl tensor 4>4, which in terms of spin- 
weight -2 spherical harmonics _2 Y* m (0, [48] reads [see 

Ref. [35] for details] 

MRtU = Y,-2Cim{t)-2Y tm {0,<p), (19) 

lm 

being R the extraction radius. In terms of the + and x 
GW polarizations we have 

^4 — — (h+ — ih x ) • ( 20 ) 

Thus, we can write 

MRCim = MR J dSl -. 2 Y; m {0,<p)(h + -ih x ). (21) 

In the adiabatic approximation (u/u 2 <Cl), we obtain 
Cim - ~m 2 u 2 hi m , (22) 


The RD modes are attached at the time when the orbital 
frequency reaches its maximum and this occurs slightly 
before the light-ring position, r ma t c h =??, Mu =?? ( rj = 
0.25) and r matc h =??, Mu =?? ( rj ~ 0.16). The matching 
is obtained by imposing the continuity of hi m (t) and all 


where 

him — ih x )im = I (^» { P) (h-\- ih x ) • 

(23) 

We compute the EOB hi m in the so-called restricted ap- 
proximation, i.e., at leading order in the PN expansion. 
They read: 


EOB 
'22 — 

-8^1 77 (Mu;) 2 / 3 e- 2iv , 

(24) 

EOB _ 
'21 - i " 

8 fn dm ... . _ iu> 

-3 V5 M’ l(M “ )e ' 

(25) 

EOB 
'32 — 

-3 ^ ( Mw ) 4/3e_2 * ¥ ’> 

(26) 

EOB 
'33 — 

- 3 \/f 

(27) 

EOB 
'44 — 

-? /y ??(1_3?7) ( Mw ) 4/3e ~ 4i< "’ 

(28) 


where 8m = m\ — m2. 

The equal-mass run lasts for ~ 14 GW cycles be- 
fore merger. This run was published originally in 
Refs. [9, 35], and further studied for data-analysis pur- 
poses in Ref. [10]. The unequal-mass runs, mi/m2.= 
2/3, 1/2, 1/4, last for 5, 5, and 9 GW cycles before 
merger. The mi/ m 2 = 2/3, 1/2 cases were published 
in References [10, 11] and the mi/ m 2 = 1/4 has recently 
been computed by the NASA-Goddard group. Adap- 
tive mesh refinement was employed for this case as in 
the previous simulations, with a finest mesh resolution of 
hf = 3M/160 used in one run and hf — M/64 used in 
a second. Adequate convergence of the Hamiltonian and 
momentum constraints were found; the numerical details 
will be reported in a future publication. In the follow- 
ing analysis of the mi /rri2 = 1/4 case, data from the 
hf = M/64 run are used. 

For different mass ratios, we show in Fig. 4 how the FF, 
obtained using white noise, depends on the parameter w. 
Based on this plot we decide to use in the rest of the 
paper w — 60. In Figs. 2 and 3, we plot the EOB and NR 
waveforms, orbital frequencies, and phase differences for 
the most accurate, long numerical simulations, notably 
the 1 : 1 and 1 : 4 cases. In Appendix B we show similar 
plots for the shorter runs 1 : 2 and 2:3. 


the higher time derivatives needed to fix the 2 N unknown 
amplitudes and phases of the N RD modes [25]. Follow- 
ing Ref. [8], we attach the fundamental mode, and two 
overtones. We find that the matching-performance does 
not improve significantly if we add more overtones. The 
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V 

\Mf / M~\ fromE 

\af / M / ] from RD 

[Mf / M~\ from fit 

[af / M /]from fit 

0.25 

0.955 

0.687 

0.958 

0.692 

0.24 

0.962 

0.671 

0.960 

0.671 

0.22 

0.967 

0.641 

0.966 

0.627 

0.16 

0.978 

0.478 

0.979 

0.482 

0.12 



0.986 

0.374 

0.08 



0.992 

0.258 

0.04 



0.997 

0.133 

0.01 



0.999 

0.034 


TABLE I: For several values of 77, we list in the second and third columns the best - values of Mf/M and af/Mf computed from 
the energy released and by extracting the fundamental QNM from - 2 C 22 , respectively. In the fourth and fifth columns we list 
the values obtained from the fits Mf/M = 1 — 0.0610 77 — 0.4339 r? 2 , and af/Mf — 3.441 77 — 2.687 77 s , respectively. 


frequency and decay time of the RD modes are computed 
using the mass and spin of the final BH, as predicted by 
the numerical simulations. In fact for non-spinning bi- 
nary systems, it is now possible to determine how the 
final-BH mass and spin depend on the mass ratio. 

In Table I we list the final BH masses and spins for 
77 = 0.25,0.24,0.22,0.16, extracted from the NR simu- 
lations. The values are compatible with Ref. [47]. We 
also computed the final BH masses and spins to lower 
values of 77 applying a two-parameter fit to the data 
77 = 0.25,0.24,0.22,0.16. For the BH mass we use the 
fit function Mf/M = 1 4- arj -f bq 2 , and get a = —0.0610 
and b = —0.4339. Quite interestingly, the value of a 
is quite close to the LSO energy for a test particle in 
Schwarzschild, i.e., i/8/9 — 1 ~ —0.0572. For the BH 
spin we employ the fit function af/Mf = cq + dr/ 2 , and 
get c = 3.441 and d = —2.687. Again, the value of c 
is quite close to the LSO angular- momentum for a test 
particle in Schwarzschild, i.e., y/l2 ~ 3.464. The ex- 


trapolation to smaller values of 77 are consistent with the 
values obtained in Ref. [45] using the EOB approach. 

To measure the differences between the NR and EOB 
waveforms we compute the fitting factor (FF), or am- 
biguity function [10, 28, 42]. The overlap (hi(t), Ji2(t)) 
between the waveforms h\{t) and /12W is defined by: 

M*).fe 2 (*)) = 4Re j ( 00 (29) 

where hi(f) is the Fourier transform of hi(t), and Sh(f) 
is the detector’s power spectral density (PSD). The FF is 
the normalized overlap between the NR waveform h NR (t) 
(target) and the EOB waveform h EOB (to, p>o) (template) 
maximized only over the initial time to and initial phase 
< po , and minimized over the initial phase of the target, 

that is 


J 


FF = min max 


(/i NR (<p ; A* ) , /i EOB (t 0 , (p 0 ; V ) ) 

\J (/i nr (<£; A*), /i nr (v?; A i ))</i EOB (i 0 , <Po; A 4 ), /i EOB (i 0 , <Po‘, A*)) 


(30) 


For the detector PSD we shall consider either white noise 
or the LIGO noise. 

In Tables II, III we list the FFs and the phase differ- 
ence (in one GW cycle) between the p4PN-EOB wave- 
form and the NR waveform, for white PSD and LIGO 
PSD, respectively. We consider the l — 2,m — 2 mode, 
but also other modes which are no-longer subdominant 
when the mass ratios decrease, as can be seen in right 
panel of Fig. 7. The FFs are rather high except for the 
l = 4 and m — 4 mode with mass-ratio 1:4. In this 
case we find that the matching procedure is not so effi- 
cient in reproducing, especially, the amplitude of the NR 
ringdown waveform. More studies extended to different 
mass ratios may shed light on this apparent anomaly. 

The hi m waveforms shown in Figs. 2, 3, 5 and 6, are 


normalized waveforms. In the right panel of Fig. 7, 
we compare the NR amplitude and EOB restricted- 
approximation amplitude for several l and m modes. 
By restricted- approximation amplitude we mean that for 
each l and m we restrict ourselves to the leading order 
term in a PN expansion. We notice, as already pointed 
out in Ref. [10], a non-negligible difference between the 
NR and EOB amplitude. We do not know what is the 
source of this difference. Higher-order PN corrections in 
the amplitude could in part nullify the difference, but 
due to the oscillatory behaviour of the higher PN correc- 
tions [8], it is hard to draw a robust conclusion. More- 
over, there could be a systematic error in the numerical 
amplitude due to extraction radius or resolution issues. 
Although, the difference in the amplitudes has negligible 
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V 

Im 


A<£gw/(27t) 

0.25 

22 

0.9907 

±0.030 

0.24 

22 

0.9881 

±0.058 

0.22 

22 

0.9878 

±0.078 

0.16 

22 j 

0.9925 

±0.035 

0.16 

33 

0.9860 

±0.055 

0.16 

44 

0.9436 

±0.065 


TABLE II: For several mass configurations, we list the FF 
obtained using white noise and maximizing only on the time 
of arrival and initial phase, and the phase difference in one 
GW cycles. 


(mi ± ra 2 ) 

lm 


A^gw/ (27t) 

(15 + 15) Mq 

22 

0.9975 

±0.043 

(50 + 50)M© 

22 

0.9817 

±0.043 

(20 + 30) M© 

22 

0.9897 

±0.065 

(20 + 40) M© 

22 

0.9889 

±0.068 

(10 + 40) M© 

22 

0.9961 

±0.035 

(10 + 40) M© 

33 

0.9911 

±0.055 

(10 + 40)M© 

44 

0.9720 

±0.075 

(20 + 80) M© 

22 

0.9965 

±0.035 

(20 + 80) M© 

33 

0.9873 

±0.055 

(20 + 80) M© 

44 

0.9548 

±0.065 


TABLE III: For several mass configurations, we list the FFs 
obtained using LIGO PSD and maximizing only on the time 
of arrival and initial phase, and the phase difference in one 
GW cycles. 


effects on the fitting- factors [ 10 ], it needs to be sorted 
out. 


V. CONCLUSIONS 

In this paper we started exploiting the flexibility of 
the EOB approach to build accurate waveforms for non- 
spinning binary BHs. We obtained phase differences be- 
tween the EOB and NR waveforms of less than 8 % in 
a GW cycle, for mass ratios rj = 0.25,0.24,0.22 and 
rj = 0.16. The p4PN EOB model can be used to predict 
waveforms for values of eta < 0.16, that can be tested 
against numerical simulations when accurate, long wave- 
forms will become available. 

The improved EOB model is obtained by simply 
adding a p4PN order constant coefficient (w — 60) in 
the EOB radial potential A(r) and Padeing the latter 
to guarantee the presence of the LSO and the light-ring. 
This modification provides a better agreement between 
the EOB and NR GW frequency during the transition 
inspiral(-plunge) to ringdown. The complete waveform 
is then built by attaching three QNMs at the peak of 
the EOB GW frequency. The QNM frequency and de- 
cay time are fixed by the final BH mass and spin, which 
can be predicted by extrapolating sparse NR results, ei- 
ther through a fit (see Table I) or within the EOB model 


itself [45]. 

We point out, as already done in Ref. [10], that there 
exists a non-negligible difference in the amplitude of the 
PN (EOB or Taylor-expanded-PN) and NR waveforms, 
whose origin has not yet been accounted for. It might be 
due to higher-order PN corrections in the amplitude or 
numerical systematic errors. 

Using white PSD, we found FFs >0.98, except for the 
l = 4 and m — 4 mode of mass ratio 1 : 4 which has FF 
= 0.94 (see Table II). At present time, we do not know 
what is the reason of this anamaly. It could be due to a 
failure of the matching procedure due to the presence of 
other QNMs. All FFs are achieved by maximizing only on 
the initial phase and time of arrival of the template (and 
by minimizing on the initial phase of the target). Since 
we do not maximize on the binary’s intrinsic parameters, 
that is the BH masses, the p4PN-EOB model can be 
used for parameter estimation. Depending on the binary 
total mass and mass ratios, FFs >0.98 are obtained with 
LIGO PSD (see Table III), except for the l = 4, m = 4 
mode, which has FF = 0.95-0.97, depending on the total 
binary mass. 

When maximizing on the binary masses, the p4PN- 
EOB waveform will have extremely high matching per- 
formances (FF >0.99) and can be used for coherent de- 
tection of non-spinning binary BHs. In the future we 
plan to extend this analysis to spinning, precessing bi- 
nary systems, using improved versions of the EOB model 
with spin [24, 27]. Once a complete analytic model for 
inspiral-merger-ringdown is worked out, we shall employ 
the detection schemes proposed in Refs. [41, 43] to deal 
with intrinsic and extrinsic parameters in the presence 
of spins, and provide accurate templates to be used by 
ground-based and space-based detectors. 
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APPENDIX A: COMPARISON OF MULTIPOLES 
WITH l ^ 2 AND m/2 

Here we discuss how the p4PN-EOB waveforms having 
l 7 ^ 2 and m /2 compare with the NR ones. We consider 
the dominant modes, that is l = 3, m = 3, l = 4, m = 4, 
l = 3 ? m = 2 and l = 2 ,m = 1 , and restrict the analysis 
to the case with mass ratio 1 : 4. 
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FIG. 5: Binary with mass ratio 1:4. In the left (right) panel we plot the hzz ( h 44) mode. We maximize only on the initial 
phase and time of arrival. 
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FIG. 6: Binary with mass ratio 1:4. In the left (right) panel we plot the h %2 (fei) mode. We maximize only on the initial 
phase and time of arrival. 


APPENDIX B: COMPARISON WITH BINARIES waveforms are much shorter, and the eccentricity during 

OF MASS RATIO 1 : 2 AND 2 : 3 the inspiral case is more pronounced. 

In this appendix we show how the p4PN-EOB model 
performs for mass ratios 1 : 2 and 2:3. In these cases the 
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FIG. 7: Binary with mass ratio 1 : 4. In the left panel we plot the dominant frequencies for several l and m modes. Generalizing 
Eq. (XXX) in Ref. [8], we define u>Dim = — (1/m) Im(C/ m /C f / m ) [11]. In the right panel we compare the NR amplitude and 
EOB restricted-approximation amplitude for several / and m modes. By restricted-approximation amplitude we mean that for 
each l and m we restrict ourselves to the leading order term. 
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FIG. 8: In the left (right) panel we plot the NR and EOB frequencies and amplitudes, and the phase difference between the 
EOB and NR h ,22 for a binary with mass ratio 1 : 2 (2 : 3). We maximize only on the initial phase and time of arrival. 
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